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Abstract 

In this correspondence, we present an algorithm for distributed sensor localization with noisy distance measure- 
ments (DILAND) that extends and makes the DLRE more robust. DLRE is a distributed sensor localization algorithm 
in R m (m > 1) introduced in [1], DILAND operates when (i) the communication among the sensors is noisy; (ii) the 
communication links in the network may fail with a non-zero probability; and (iii) the measurements performed to 
compute distances among the sensors are corrupted with noise. The sensors (which do not know their locations) lie 
in the convex hull of at least m + 1 anchors (nodes that know their own locations.) Under minimal assumptions on 
the connectivity and triangulation of each sensor in the network, this correspondence shows that, under the broad 
random phenomena described above, DILAND converges almost surely (a.s.) to the exact sensor locations. 
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I. Introduction 

Localization is an important problem in sensor networks, not only on its own right, but often as the first step 
toward solving more complicated and diverse network tasks, which may include environment monitoring, intrusion 
detection, and routing in geographically distributed communication networks. The problem we consider is when a 
large number of sensors do not know their locations, only a very few of them know their own. In [1], we presented 
a distributed sensor localization (DILOC) algorithm in K m (m > 1), when we can divide the ./V nodes in the 
sensor network into these two sets: the set (tofti anchors where n > m + 1 and the set of M sensors, with 
typically N » n. The n anchors are the nodes that know their exact locations, whereas the M sensors are the 
nodes that do not know their locations 1 . We assume that the sensors lie in the convex hull of the anchors, i.e., 
C(O) C C(k), where C(-) denotes the convex hull 2 . To each sensor I in the network, we associate a triangulation 
set 3 , O;, which is a set of m + 1 neighboring nodes such that sensor I lies in their convex hull, i.e., I e C(Oz). In 
DILOC, each sensor, I, updates its location estimate as a linear convex combination of the estimates of the nodes 
in its triangulation set, 0;, where the coefficients of the linear combination are the barycentric coordinates. Under 
minimal assumptions on network connectivity, and that the sensor I knows the precise distances in the set 

2} = {Z}U6 ; , (1) 

reference [1] shows that DILOC converges to the exact sensor locations. An interesting contribution of DILOC 
is that it reduces the centralized nonlinear problem of localization to a linear distributed iterative algorithm under 
broad assumptions, which can be implemented through local inter-sensor communication in real-time. 

Reference [1] extends DILOC and presents a distributed localization algorithm in random environments (DLRE). 
DLRE is a stochastic approximation version of DILOC where the DILOC iterations are weighted with a decreasing 
weight sequence that follows a persistence condition. DLRE operates under the following random phenomena: (B.l) 
the communication among the sensors and the nodes in their triangulation set is noisy, i.e., at the t-th iteration, 
sensor I receives only corrupted versions, yj n (t), 1 < j < m, of the m components of the neighboring node n's 
state, i.e., a;£(t), given by 

yin(t) = <(t) + vl(t), nee,, (2) 

where < vj n (t) > is a family of independent zero-mean random variables with finite second moments; (B.2) each 

^ J l,n,j,t 

inter-node communication link is modeled by a binary random variable, ei n (t),l G fi, n £ ®i, which is 1 (active 
link) with probability qi n , and (link failure) with probability 1 — q n i; and (B.3) the local barycentric coordinates, 
computed at iteration t from the current noisy distance measurements, can be represented as a perturbation of 
the exact barycentric coordinates. Reference [1] shows that, under (B.1)-(B.3) and if the perturbation in (B.3) is 
unbiased, DLRE converges to the sensor exact locations; however, if the perturbation in (B.3) is biased, DLRE 
converges with a steady-state error (bias). 

In this correspondence, we modify DLRE to present the algorithm DILAND (distributed sensor localization with 
noisy distance measurements) and show that it converges a.s. to the exact sensor locations under much broader 

'in the sequel, we always use this disambiguation for sensors and anchors. When the statement is true for both sensors and anchors, we use 
the term node. 

2 The minimal number of anchors required for a non-trivial convex hull in ra-dimensional (mD) space is m + 1 that is a triangle in 2D space. 
We may have more than m + 1 anchors forming the boundary of a polygon for less stringent requirements on sensor placement, see [2] for 
details. 

3 In [1], we give a convex hull inclusion test to verify a triangulation set, and we relate the communication radii and the density of deployment 
to guarantee triangulation with a high probability. We also study the probability of a successful triangulation at each sensor. For this study, we 
assume a randomly deployed sensor network where the number of sensors in any given area follows a Poisson distribution. 
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distance measurement noise assumptions. In DILAND, we replace (B.3) above with the following weaker condition: 
(B.3) at every iteration t, we assume there exist computationally efficient estimates of the required inter-sensor 
distances based on all distance measurements till time t such that these estimates are consistent, i.e., they converge 
a.s. to the exact distances as t — > oo. The state update in DILAND uses these estimates to compute the local 
barycentric coordinates, whereas the update in DLRE uses only the current distance measurements. The consistency 
assumption on the estimates of the inter-sensor distances is quite weak and, as will be shown, is applicable under 
practical schemes of estimating inter-sensor distances through: (i) received signal strength (RSS) and (ii) time-of- 
arrival (TOA) (see [3].) We emphasize that DILAND does not require spatial or temporal distributional assumptions 
on either the communication or the distance measurement noises, except for finiteness of the second order moment, 
see (B.l). 

Because of (B.3), DLRE [1] converges to the exact sensor locations with a steady state error when the resulting 
perturbation of the barycentric coordinates is biased. In contrast, under the new assumption (B.3), DILAND 
converges to the exact sensor locations regardless of the bias introduced in the system matrix (at each iteration) 
due to noisy distance measurements. This new setup leads to a behavior and analysis for DILAND that is different 
from DLRE's in [1]. This is because using distance estimates based on the entire past leads to an inherent strong 
statistical dependence in the iterative scheme. This dependence makes the analysis of DILAND different from 
standard stochastic approximation arguments, which were used in [1] to prove the convergence properties of DLRE. 
Furthermore (as we will show in Section III-C), if we do not have link failures and communication noise, the weight 
sequence, a(t), in DILAND does not require the square summability condition required by DLRE. Hence, DILAND 
can be designed to converge faster than DLRE by choosing the DILAND weight sequence to sum faster to infinity 
than the DLRE weight sequence. 

We describe the rest of the paper. Section II briefly recapitulates our prior work on DILOC and DLRE. We then 
study the distance estimates and present the main result, DILAND, of this correspondence in Section III. Section IV 
presents simulations, and Section V concludes the paper. 

II. Prior Work 

We briefly recapitulate the setting of the distributed localization problem, details and discussions are in [1]. 

A. Distributed Localization Algorithm (DILOC) 

Let the row vector xj(i) G R lxm ,Z G 17, be the m-dimensional state that represents the estimated coordinates 
of sensor I at time t. Similarly, let the row vector £ K lxm , k G k, be the m-dimensional state of the location 
coordinates of anchor k. DILOC updates at time t are the following: 

u fc (t + l) = u fc (i), k£n, (3) 

neQn&i kennQi 

The symbol 0; is the triangulation set of sensor I, and pi n and bu~ are the sensor-sensor and the sensor-anchor 
barycentric coordinates, [4], respectively; these are computed using the inter-node distances in T>i and the Cayley- 
Menger determinants [5] (see [1].) Let the set of inter-sensor distances required to compute all the barycentric 
coordinates be (* represents the exact value) 

d* = {d* kn \k,n£V l ,l£n}. (5) 
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The barycentric coordinates lead to the system matrices P(d*) = {pi n } and B(d*) = {b tk } of DILOC (we 
use (d*) to show an implicit dependence on the inter-sensor distances.) We define the matrix X(£) collecting the 
row vector states x/(£) for all sensors I and with column vectors x J '(t); similarly for U. 

M x m : X(t) = [xi (tf , . . . , x M (t) T ] T = [x 1 (*)■■■ x 3 ' (*)■■■ x m (t)] (6) 
(m + l)xm:U = [uf , . . . , T = [u 1 ^) • • • u'(i) • • • u TO (t)] . (7) 

In (6) and (7), sub and superindices indicate row and column vectors of the corresponding matrices, i.e., x J (i) is 
the column vector that collects component j of the row vector state x;(t) of all sensors I; similarly for u j (t). 
We recall from [1] the following assumptions. 

Structural assumptions: (Al) There are at least m + 1 anchors, i.e., \k\ = n = m + 1, that do not lie on 
a hyperplane in R m ; (A2) The M sensors lie inside the convex hull of the anchors, i.e., C(f2) C C(k), where 
C(-) denotes the convex hull; (A3) There exists a triangulation set 4 , 8j C 6, V I G il, with |9;| = m + 1, such 
that I e C(0/); (A4) The sensor I is assumed to have a communication link, n —> I, to each n £ 0/ and the 
inter-node distances among all the sensors in the set T>i are known at sensor I; and (A5) Each anchor, k € n, has 
a communication link to at least one sensor in fl 

Theorem 1 (Theorem 1, [1]) Under (A1)-(A5), for DILOC (4), lim^oo X(t + 1) = (I - P(d*))- 1 B(d*)U, i.e., 
the states converge to the exact sensor locations. 



B. Distributed Localization in Random Environments (DLRE) 



We start by contrasting assumptions (B.3) and (B.3) (introduced in Section I). 

(B.3) Small perturbation of system matrices: Recall from Section I that at each sensor I the distances required to 
compute the barycentric coordinates are the inter-node distances in the set T>i. In reality, the sensors do not know the 
precise distances, d^At.y but estimate the distances from RSS or TOA measurements at time t. When we 

have noisy distance measurements, we iterate with system matrices P(d t ) and B(d t ), not with P(d*) and B(d*), 
where d t is the set like d* in (5) that collects all required inter-node distance estimates dk n (t), k,n € T>i, I G O. 
Assuming the distance measurements are statistically independent over time, P(d t ) and B(d t ) can be written as 

P (d t ) = P(d*) + S P + S P (t) 4 {p ln (t)}, B (d t ) =B(d*) + S B + S B (t) = {bin(t)}, (8) 

where Sp and Sb are mean measurement errors, and {Sp(i)} t > and {SB(i)}t>o are independent sequence 
of random matrices with zero-mean and finite second moments. In particular, even if the distance estimates are 
unbiased, the computed P(d t ) and B(d t ) may have non-zero biases, Sp and Sb, respectively. The following is 
the small bias assumption, we made in [1]: p(P(d*) + Sp) < 1. The DLRE algorithm and its convergence are 
summarized in the following theorem. 

Theorem 2 (Theorem 3, [1 ]) Under the noise model (B.1)-(B.3), the DLRE algorithm given by 

x /( y + i)= (i-M(/))x,c/) + M(n J2 einmn[t) (x n (t) + v, w (t))+ £ e '*(*)M*) (Ufe + v<fe(f)) 

nenne, qin fce K ne, qik 

(9) 



In [1], we give a convex hull inclusion test to identify such triangulation set at each sensor /. 
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for I G 0, with a(t) satisfying the persistence condition a(t) > 0, ^ t a(t) = oo, J2t a ' '(*) < °°> converges almost 
surely to lim^ X(i + 1) = (I - P - S P )- 1 (B + S B )U(0). 

The DLRE converges to the exact sensor locations for unbiased random system matrices, i.e., Sp = Sb = 0, as 
established in Theorem 2. As pointed out earlier, even if the distance estimates are unbiased, the system matrices 
computed from them may be biased. In such a situation, the DLRE leads to a nonzero steady state error (bias). 

III. Distributed Sensor Localization with Noisy Distance measurements 

As mentioned, since DLRE at time t uses only the current RSS or TOA measurements to compute distance 
estimates, the resulting system matrices have an error bias, i.e., Sp ^ and Sb ^ 0. We use the information from 
past distance measurements to compute the system matrices at time t; these become a function of the entire past 
measurements, {RSS s } s < t or {TOA s } s < t - The DILAND algorithm efficiently utilizes the past information and, 
as will be shown, leads to a.s. convergence to the exact sensor locations under practical distance measurement 
schemes in sensor networks. To this aim, we review typical models of distance measurements in wireless settings 
in Section III-A and introduce the DILAND algorithm in Section III-B. Section III-C discusses DILAND. 



A. Models for distance measurements 

We explore two standard sensor networks distance measurements: Received Signal Strength (RSS) and Time-of- 
Arrival (TOA). We borrow experimental and theoretical results from [3]. 

1) Received signal strength (RSS): In wireless, the signal power decays with a path-loss exponent, n p , which 
depends on the environment. If sensor a sends a packet to sensor b, then RSS a & is the power of the signal received 
by sensor b and the maximum likelihood estimator, d ab , of the distance, d ab , between sensors a and b is [3] 

~ n -Rss ab 

<U = A 10 10 " P , (10) 

where n is the received power at a short reference distance A . For this estimate, 

E[d ab \=Cd ab , (11) 

where C is a multiplicative bias factor. Based on calibration experiments and a priori knowledge of the environment, 
we can obtain precise estimates of C; for typical channels, C « 1.2 [6] and hence scaling (10) by C gives us an 
unbiased estimate. If the estimate of C is not acceptable, we can employ the following scheme. 

DILAND is iterative and data is exchanged at each iteration t. We then have the measurements on RSS and C 
at each iteration. Hence, the distance estimate we employ is 

d ab (t) = (12) 

where C(t) is the estimate of C (for details on this estimate, see [3] and references therein) at the t-\h iteration 
of DILAND. Assuming that d ab (t) and C(t) are statistically independent (which is reasonable if we use different 
measurements for these estimates and assume that the measurement noise is independent over time), we have 



E d ab (t) 



dab- (13) 



Since at time t, we have knowledge of {d a b{s)}o<s<t, we can use the following sequence, {d a b(t)}t>o, of 
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distance estimates to compute the barycentric coordinates at time t: 

d ab (t) = \Yj = t -^~d ab (t - 1) + jd ab (t), d ab (0) = d ab (0). (14) 

s<t 

Then, from (12)-(13) and the strong law of large numbers, we have d ab (t) — > d ab a.s. as i — > oo. 

2) Time-of-arrival (TOA): Time-of-arrival is also used in wireless settings to estimate distances. TOA is the time 
for a signal to propagate from sensor a to sensor b. To get the distance, TOA is multiplied by v p , the propagation 
velocity. Over short ranges, the measured time delay, T ab , can be modeled as a Gaussian distribution 5 [3] with 
mean d ab /v p + /iy and variance a\. The distance estimate is given by d ab = (T ab — /j,t)v p . Based on calibration 
experiments and a priori knowledge of the environment, we can obtain precise estimates of the bias ht\ wideband 
DS-SS measurements [7] have shown [It = 10.9 ns. 

Since DILAND is iterative, we can make the required measurements at each iteration t and compute an estimate, 
juT(i), of the bias, ht (for details on this computation, see [7]). Then, using the same procedure described for RSS 
measurements, we can obtain a sequence, {d ab (t)} t >o, of distance estimates such that d ab (t) — > d ab a.s. as t — > oo. 

In both cases, we note that, if {Z(t)} t >a is a sequence of distance measurements, where Z = RSS or Z = TOA, 
collected over time, then there exist estimates d t with the property: d t — > d* a.s. as t — > oo. In other words, by 
a computationally efficient process (e.g., simple averaging) using past distance information, we can estimate the 
required inter-sensor distances to arbitrary precision as t — > oo. This leads to the following natural assumption. 

(B.3) Noisy distance measurements: Let {Z(t)} t >o be any sequence of inter-node distance measurements 
collected over time. Then, there exists a sequence of estimates {d t } t >o such that, for all t, d t can be computed 
ejficientiy from {X(s)} s < t and we have 

'lim d t = d*l = 1 (15) 



We now present the algorithm DILAND under assumption (B.3). 
B. Algorithm 

For clarity of presentation, we analyze DILAND in the context of noisy distance measurements only and assume 
that the inter-sensor communication is perfect (i.e., no link failures nor communication noise 6 .) Let P(dt) = {Pi n (t)} 
and B(d t ) = \bi k (t)} be the matrices of barycentric coordinates computed at time t from the distance estimate d t . 
The DILAND algorithm updates the j-th component of the state of all sensors, i.e., updates x J (t), the j-th component 
of the location estimates of all sensors, as follows: 

x*(i + l) = {l-a(t))x?(t)+a(t) [P(d t )x^(i)+B(d t )u J ] , 1 < j < m, (16) 

where the weight sequence a(t) satisfies a(t) > 0, lim^oo a(t) = 0, and ^ t a(t) = oo. In particular, here we 
consider the following choice: for a > and < 5 < 1, 

The update (16) is followed by the distance update (14). The following result gives the convergence properties of 
DILAND. The proof is provided in Appendix I. 



5 Although, as noted before, DILAND does not require any distributional assumption. 
6 In Subsection III-C, we discuss the effect of link failures and communication noise. 
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Theorem 3 Assume (A.1)-(A.5) and (B.3). Let {x J (i), 1 < j < m} t >o be the j-th coordinate of the state sequence 
generated by DILAND (16). Then x J (i) converges a.s. to the exact j-th coordinate location as t — > oo, i.e., 



lim *>{t) = (I - P(d*)) _1 B(d*)u J ', Vj = 1, . . . 



rri 



1, 



(18) 



which are the exact sensor locations as established in Theorem 1. 



C. DILAND.Discussions 

We discuss the consequences of Theorem 3. Unlike DLRE, the proof of Theorem 3 does not fall under the 
purview of standard stochastic approximation. This is because the system matrices P(d t ), B(d t ) at any time t are 
a function of past distance measurements, making the sequences of system matrices a strongly dependent sequence. 
On the contrary, DLRE assumes that the sequences of system matrices are independent over time. Thus, DLRE 
can be analyzed in the framework of standard stochastic approximation (see, for example, [8]), where it is assumed 
that the random perturbations are independent over time (or more generally martingale difference sequences.) In 
Appendix I, we provide a detailed proof of Theorem 3, which also provides a framework for analyzing such 
stochastic iterative schemes with non Markovian perturbations. 

For clarity of presentation, we ignored the effect of link failures and additive communication noise in the analysis 
of DILAND. In the presence of such effects (i.e., assumptions (B.1)-(B.2)), DILAND can be modified accordingly 
and the corresponding state update equation will be given by: for 1 < j < m 

x»'(t + 1) = (1 - a(t))^(t) + a(t) [E P (d t ) (x»'(t) + v j \t)) + E B (d t ) (u j + v j (i))] , 1 < j < m, (19) 

where v J (t) is the measurement noise under (B.l), is the Hadamard or pointwise product of two matrices, and 
= | | models the link failures as per (B.2), see Section I. Here we use the notation of Sections I and II in 
the context of random environments. However, for (19), the weight sequence {a(t)} needs to satisfy the additional 
square summability condition, i.e., J2t>o a2 (t) < °°- With this modification, the results in Theorem 3 will continue 
to hold. The proof will incur an additional step, i.e, after establishing Theorem 3 with no link failures and noise, 
a comparison argument will yield the results (see [9] for such an argument in the context of distributed estimation 
in sensor networks.) In fact, considering a parallel update scheme with state sequence {x(t)} given by 

x*'(t + 1) = (1 - a(t))5? (t) + a(t) [P (d t ) 5? (t) + B (d t ) u j ] , 1 < j < m, (20) 

we note that the sequence {x 3 (t)} converges a.s. to the exact sensor locations by Theorem 3. Subtracting this 
fictitious update scheme from the modified DILAND iterations in eqn. (19), we note that the residual eP (t) = 
x- 7 (t) — x 3 (i) evolves according to a stable system driven by martingale difference noise (here we note that the 
effects of link failures and additive channel noise are considered to be temporally independent.) The a.s. convergence 
of the sequence {e J (t)} to then follows by standard stochastic approximation arguments (see, for example, [8].) 

In the absence of link failures and additive communication noise (possible with a TCP-type protocol and digital 
communication,) the DILAND needs the weight sequence a(t) to satisfy the weaker condition on the weights below 
(16) rather than J2t>o a2 (t) < 00 • ^ n tne contrary, the DLRE (even in the absence of link failures and additive 
communication noise) requires the additional square summability condition on the weight sequence, because of 
the presence of independent perturbations {Sp(t)} t >o and {SB(i)}t>o with non-zero variance. Thus, DILAND 
can be designed to converge faster (see also Section IV for a numerical example in this regard) than DLRE by 
choosing a weight sequence that sums to infinity faster than can be achieved in the DLRE. This is because the 
convergence rate of such iterative algorithms, essentially, depends on the rate at which the iteration weights sum to 
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infinity. In particular, the square summability assumption on the weight sequence as required by DLRE ([1]) limits 
the convergence rate to the order of ([9]), whereas the DILAND does not require square summability of the 
weight sequence and much better convergence rates can be obtained by making {a(t)} sum to oo faster. However, 
the requirement that a(t) — > as t — > limits the convergence rate achievable by the DILAND, which is always 
sub-exponential as exponential convergence rates require the a(t) to be bounded away from zero in the limit. 7 

With respect to triangulation, the variance of the distance measurements may be so large that the convex hull 
inclusion test provided in [1] may lead to incorrect triangulations in the DLRE. However, DILAND refines the 
inter-node distance estimates at every time-step, so eventually these estimates get more accurate and one will end 
up with the correct triangulation, possibly, with a different set at each iteration. Interested readers can refer to [2] 
for a localization algorithm when the triangulation set, Qi, is a function of t, i.e., 0j(t). 

Remark. DILAND estimates and updates the inter-node distances at each time step. An alternative is to estimate 
well the inter-node distances by initially collecting and averaging several measurements and then to run DLRE. 
There are clear tradeoffs. The second scheme introduces a set-up phase, akin to a "training phase," implying a delay 
in getting distance information, it is sensitive to variations in sensor locations, and it will show a remnant bias 
because it stops the estimation of the inter-node distances after a finite number of iterations, which means there will 
be a residual error in the inter-sensor distances. This residual error induces a bias in the distance estimates and so 
in the barycentric coordinates, which means a residual error in the sensor coordinates, no matter how long DLRE is 
run. DILAND does not suffer from these limitations, it has no initial delay, it is robust to variations in the sensors 
geometries since it can sustain these to a certain degree given its adaptive nature, and it continues reducing the 
asymptotic error since the inter-sensor distances are continuously updated. On the other hand, DILAND's iterations 
are more onerous since it requires continuous updating of the inter-node distances and recomputing the barycentric 
coordinates. Very similar tradeoffs arise an are studied in the context of standard consensus in [10] 

IV. Simulations 

In this section, we compare DLRE [1] with DILAND. In the first simulation, we do not include link failures and 
communication noise and focus only on distance noise. This study presents the advantage of choosing the relaxed 
weight sequence, a DILAND (t) below (16), over a DLRE in Theorem 2. We assume that the distance, dij, between any 
two sensors I and j is corrupted by additive Gaussian noise, i.e., 

d lj (t)=dt j +w lj (t), (21) 

where tuy ~ Af (o, We choose the variance of the noise in distance measurement as 10% of the actual 
distance. We adopt the estimator (14), which is an estimator at time t given by the average of the past observations. 
Hence, (14) becomes our distance update. It is straightforward to show that dij(t) — > d\-. We use the following 
weight sequences for DLRE and DILAND compared in Fig. 1 (center): 

Q! DLRe(^) = Q 55 ) C^DILAND^) = q 25 • (22) 

These satisfy the persistence condition in Theorem 2 and below (16), respectively. In particular, as noted before, 
«diland does not have to satisfy J2t>a "dilandW < °°> which leads to faster convergence (see Fig. 1 (right.)) 

We simulate an N = 50 node network inm = 2-D Euclidean space. We have m + 1 = 3 anchors (with known 
locations) and M = 47 sensors (with unknown locations). This network (and appropriate triangulations) is shown 

7 The assumption a(t) — > cannot be relaxed in general. This is beacuse, although the distance estimator is assumed to be consistent, at 
every time t, there is an estimation error, and these errors may accumulate in the long run if a constant a(t) was used. 
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Fig. 1. Distance noise only: Left: TV = 50 node network, m = 2-D space, m + 1 = 3, anchors and M = 47 sensors. Center: Decreasing 
weight sequences, ct(t), for DLRE and DILAND. Right: Normalized mean squared error for DLRE and DILAND. 




— DILAND 
— DLPL ;i 



1000 1500 2000 2500 
Iterations, t 



Fig. 2. Distance noise, communication noise, and link failure: Left: N = 50 node network, m = 2-D space, m + 1 = 3 anchors, and M = 47 
sensors. Center: Same weight sequence, a(t), for DLRE and DILAND. Right: Normalized mean squared error for DLRE and DILAND. 



in Fig. 1 (left), where Vs represent the anchors and os represent the sensors. Fig. 1 (center) shows the weight 
sequences chosen for DLRE and DILAND, whereas Fig. 1 (right) shows the normalized mean squared error 

MSE t = M > (23) 

jjYZiY^i (*?(o)-»T) 

In our second simulation, we keep the experimental setup, but include a zero-mean Gaussian random variable with 
unit variance as communication noise. We further assume that the communication links are active 90% of the time. 
With this additional randomness in the localization algorithm, we cannot use the relaxed DILAND weights, but use 
the same weights given in Theorem 2 for both DLRE and DILAND, see Section III-C on a discussion on this choice. 
The network (along with appropriate triangulations) is shown in Fig. 2 (left.) Fig. 2 (center) shows the weights 
sequences chosen for DLRE and DILAND, whereas Fig. 2 (right) shows the normalized mean squared error. In both 
studies, we note that DLRE converges with a steady state error as is emphasized earlier, whereas, DILAND converges 
almost surely to the exact sensor locations. Furthermore, the simulations confirm that DILAND converges faster 
(faster convergence in the presence of all random phenomena is a consequence of refined distance measurements) 
because the square summability condition on a DLRE (t) as required by DLRE is not required by a D1LAND (t). 
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V. Conclusions 

In this correspondence, we present a distributed localization algorithm, DILAND, that converges a.s. to the 
exact sensor locations when we have (i) communication noise, (ii) random link failures, and (iii) noisy distance 
measurements. We build on our earlier work, DLRE [1], that assumes that the noise on distance measurements 
results into a small signal perturbation of the system matrices; this perturbation is biased, in general. Due to this 
bias, there is a non-zero steady error in the location estimates of DLRE. With the additional assumption that the 
distance estimates at each iteration of DILAND converge a.s. to the exact distances, we show that the steady state 
error in the location estimates of DILAND is zero. We further show that if there is no communication noise and 
link failures, possible in some communication environments, the convergence rate of DILAND is faster than the 
convergence rate of DLRE. We provide simulations to assert the analytical results. 

Appendix I 
Convergence of DILAND 

We first summarize a relevant result that is needed for proving Theorem 3. The proof of this Lemma is in [9]. 



Lemma 1 (Lemma 18, [9]) Let the sequences {ri(t)} t > and {r2(i)} t >o be given by 

ffll ,, N 0,2 



ri(t) 



r 2 (t) = 



where oi, 02, £2 > and < Si < 1. Then, if 5i = 82, there exists K > 0, such that, for non-negative integers, s < t, 



II (l-n(0) 



=fc+l 



r 2 (k) < K 



(24) 



(25) 



Moreover, the constant K can be chosen independently of s,t. Also, if 61 < <5 2 , then, for arbitrary fixed s, 



lim y 



II (l-n(0) 



=fc+l 



r 2 (k) = 



We now prove the convergence of DILAND given by (16). Unless otherwise noted, the norm 
Euclidean 2-norm. The following lemma shows that the DILAND iterations are bounded for all t. 



(26) 

refers to the standard 



Lemma 2 Consider the sequence of iterations in (16). We have 

P sup ||x i (i)|| < 00, l<j<m 
_t>o II II 

In other words, the sequence {x J (i)} t > remains bounded a.s. for all j. 
Proof: We rewrite (16) as 

x 3 (t + l) = (l-a(t))x j (t) + a(t) [(P(d t )+P(d*)-P(d*))x 3 (t) + B(di)u J ] , 

= [(1 - a(t)I) + a(i)P(d*)] x J (t) + a(t) [P(d t ) - P(d*)l x 3 (t) + a(i)B(d t )u 



(27) 



(28) 
(29) 



Since p(P(d*)) < 1 (recall Theorem 1), it follows from the properties of linear operators on Banach spaces that there exists 
a norm 



\ p such that the corresponding induced norm of the linear operator P(d*), satisfies 

l|P(d*)y|| P 



|P(d* 



max ■ 



llyllp 



< 1 



(30) 



Moreover, such a norm can be chosen to be equivalent to the Euclidean norm ||-|j (see, for example, [11]), i.e., there exist 



constants ci,c 2 > 0, such that, ci ||-|| < 



< C 2 



In other words, 



and 



generate the same topology on the 
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Euclidean space. Since a(i) — + 0, we can choose to sufficiently large, such that, a(t) < 1, Vt > to- From the above 
construction, we have for, t > to, 



(1 - a(t)I) + a(t)P(d*)|| P < 1 - a(t) + a(t) ||P(d*)|| P = 1 - A*a(t) 



(31) 



where < A* < 1. From (B.3), we recall d t — > d* a.s. and from P(-), B(-) being continuous functions of d t , we have 

P(d t ) P(d*), a.s. B(d t ) B(d*), a.s. (32) 

Now fix a sample path ui. There exists ti(a>) sufficiently large, such that, for < e < A*, we have, if t > ti(u>) 



|P(d t )-P(d*)|L< e 



(33) 



Also, the a.s. convergence of B(d t ), implies there exists Ai(w) such that ||S(d t )|| < Ai(u>). Then, for t > max(£ ,ii(w)) 



x 3 (i + l) < (1-A*) x J (i) + ea(t) x J (t) + a(t)Ai(u>) \\u J \\ 
I Hp II Hp H Hp H Hp 

= -e)a(t))\\^(t)\\ +a(t)A 1 (^)||u 3 || 

Hp Hp 



(34) 



Let t2{u>) = max(to, ti(u>)), ai — A* — e, and 02(0;) = Ai(cj) lu^ II . Continuing the above recursion, we have for t > t2{u>) 



W(t)\\ < 
Hp 



J] (l-aia(«)) ||x(t 2 H)|| p + £ 

vfc=t 2 (w) / fc=t 2 (w) 



]{ (1 - oia(/)) a 2 a(k) 



(=fe+i 



t-i 



< ||x(t 2 H)|| P + 

k—t2 (ijJ 



]^[ (l-oia(I)))o2(w)a(*) 

l=k+l 



The second term falls under the purview of Lemma 1 in Appendix I with 5i = 82 = 8, and we have 

sup ||x J (t)j| <K(u>) 



(35) 



(36) 



for some K{uS) > 0. Since the above analysis holds a.s., we have P [sup t>0 ||x J (t)|| < 00] =1. The lemma then follows 
from the fact that the norms ||-|| and ||-|| p are equivalent. In particular, boundedness in ||-|| is equivalent to boundedness in ||-|| p . 

■ 

We now present the proof of Theorem 3. 
Proof: [Proof of Theorem 3] We use a comparison argument. To this end, consider the idealized update 

Z j (t + 1) = (1 -a(t))5c>(t) +a{t) [P(d*)x(t) +B(d> J ] (37) 

It follows from Theorem 1 that, for j — 1, . . . , m, 

lim x'(t) = (I - P(d*))" 1 B(d*)u J . (38) 

Define the sequence {e j (t) = x J (t) — x J (t)} t >o. We then have 

e J (t + 1) = (1 - a(t)) e 3 (t) + a(t)P(d*)e 3 (t) + a(t) (P(d t ) - P(d*)) x J (i) + a(t) (B(d t ) - P(d*)) u J . (39) 

Recall the norm ||-|| p introduced in the proof of Lemma 2. From (31), we have for t > to, 

||(1 - a(i)I) + a(t)P(d*)\\ p < 1 - \*a{t) (40) 

where < A* < 1. Now, fix a sample path (lo). Consider r > 0. By (32), there exists sufficiently large £3(0^7-) such that 



Vt>t 3 (w,T): ||B(d t )-B(d*)|L<T, 



|P(d«)-P(d*)|L<T 



(41) 
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Define £4(1^, r) = max t0it3 ( UiT ). We then have for t > ti(uj, r) 
||e J '(t + l)|| 



< 


11(1- 


-a(t)7) + a (t)P(d')|| P e J (t) \\ p 




+ 


a(t) 


||P(d t )-P(d*)|| p x J (t) || p + a(*) 


|B(d t )-B(d*)|| p 


< 


(1- 


A*a(t)) e J '(t) 1 +Ta(t)/f(w) + r 


\\p 


< 


(1- 


A*a(t)) \\e j (t) \\ p +ra(t) (k{u) + \ 


U l) 



(42) 



(43) 



where K{w) is defined in (36). Continuing the above recursion, we have for t > £4(^,7-) 

H (l-A'a(O)) (*(a>) + |u'|Ja(*) 



H*)||p^ II (1-A'a(*))]||e''(t4(w,r))|| p + T £ 

£4(^,-7-) / fc— £4(cj,tJ 



i=fe+l 



Using the inequality 1 — X*a(t) < e x for sufficiently small X*a(t), we have 

TT (1- X*a(t)) I ||e 3 (t 4 (cj,r))|| < lim e^* E *='4(-.') a(t) ) 

J- J- /II Hp t— *oo 



lim 

t — >oo 



1 k — £4 (tjj ,t) 

The last step follows from the fact that YlkLt (u r) a W = 00 ■ From Lemma 1, we have 



|e J (t 4 (o;,T))|| =0 



k = t4(uj,r) 



TJ (l-A*a(J))J (*(<") + 



(=fe+i 



a(jfc) 



< c 4 (cj) 



(44) 



(45) 



(46) 



Note, in particular, 04(0;) is independent of r. We then have from (44) 



limsup||e J (t)| < limsup ] [ (1 - X*a(t)) \ \\e J (t 4 (uj, r))| 



+ r lim sup 



fc=£ 4 (uj ,t) 



(l-A'a(O)J (5H + ||u J ||Ja(fc) 



< rc4(a»). 



Since (47) holds for arbitrary r > 0, we have lim^oo ||e J (i)|| = 0. The theorem then follows from the fact that this 
convergence to zero holds for u> a.s and from the equivalence of the norms ||-|| and ||-|| p . ■ 
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